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ABSTRACT 

We investigate the accretion of angular momentum onto a protoplanet system 
using three-dimensional hydrodynamical simulations. We consider a local region 
around a protoplanet in a protoplanetary disk with sufficient spatial resolution. 
We describe the structure of the gas flow onto and around the protoplanet in 
detail. We find that the gas flows onto the protoplanet system in the vertical 
direction crossing the shock front near the Hill radius of the protoplanet, which 
is qualitatively different from the picture established by two-dimensional simu- 
lations. The specific angular momentum of the gas accreted by the protoplanet 
system increases with the protoplanet mass. At Jovian orbit, when the proto- 
planet mass Mp is Mp < IMj, where Mj is Jovian mass, the specific angular 
momentum increases as j oc Mp. On the other hand, it increases as j oc 
when the protoplanet mass is Mp > IMj. The stronger dependence of the specific 
angular momentum on the protoplanet mass for Mp < IMj is due to thermal 
pressure of the gas. The estimated total angular momentum of a system of a gas 
giant planet and a circumplanetary disk is two-orders of magnitude larger than 
those of the present gas giant planets in the solar system. A large fraction of 
the total angular momentum contributes to the formation of the circumplanetary 
disk. We also discuss the satellite formation from the circumplanetary disk. 

Subject headings: accretion, accretion disks — hydrodynamics — planetary sys- 
tems — planets and satellites: formation — solar system: formation 
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INTRODUCTION 



Until now, more than 200 extrasolar planets (or exoplanets) have been detected mainly 
by measuring the radial motion of their parent star along the line of sight. Almost all 
exoplanets observed by this method are giant planets, like Jupiter and Saturn in our solar 
system, because massive planets are preferentially observed. Although these planets are 
supposed to be formed in the disk surrounding the central star (i.e., circumstellar disk or 
protoplanetary dis k), their formation pro cess has not been fully understood yet. In the core 
accretion scenario (IHayashi et al. Ill985l ). a solid core or protoplanet with ~ 10 M®, where 

is the Earth mass, captures a massive gas envelope from the protoplanetary disk by 
self-gravity to become a gas giant planet. 

The evolution of the gaseous protoplanet has b een studied w i th the approximation of 
spherical symmetry including radiative transfe r 



Pollack et al. 


1996: 


Ikoma et al. 


2000 


)■ 



(e.g., 


Mizuno 


1980; 


Bodenheimer & Pollack 


Ikoma et al. 


(200C 


) showed that rapid gas 



accretion is triggered when the solid core mass exceeds ~ 5 — 20 M®, and the protoplanet 
quickly increases their mass by gas accretion. However, the angular momentum of the accret- 
ing gas was ignored in these studies, because they assumed spherical symmetry. Since the gas 
accretes onto the solid core with a certain amount of the angular momentum, a circumplan- 
etary disk forms around the protoplanet, analogously to the formation of a protoplanetary 
disk around a protostar. The difference in the disk formation between the protostar and 
protoplanet is the region from which the central object acquires the angular momentum. 
The protostar acquires the angular momentum from a parent cloud, while the protoplanet 
acquires it from the shearing motion in the protoplanetary (circumstellar) disk. In addition, 
the gravitational sphere of the protostar spreads almost infinitely, while the gravitational 
sphere (i.e., the Hill sphere) of the protoplanet is limited in the region around the proto- 
planet because the gravity of the central star exceeds that of the protoplanet outside the 
Hill sphere. 

The numerical simulations are useful to investigate the gas accretion onto a proto- 
planet and its circumplanetary disk (hereafter we just call them as a protoplanet sys- 
tem). iKorycansky et al. I (Il99ll ) studied giant planet formation using one-dimensional quasi- 
spherical approximation with angular momentum transfer. They showed that as the proto- 
planet contracts, outer layers of the envelope containing sufficient specific angular momentum 
remain in bound orbit, and form a circumplanetary disk. However, owing to the spherical 
symmetry, the accretion fiow from the protoplanetary disk to the protoplanet could not be 
investigated in their study. To study the accretion fiow onto a protoplanet and the acquisi- 
tion process o f the a ngular momentum in detail, a multidimensional simulation is necessary. 
Sekiya et al. I (119871 ) investigated the gas fiow around a protoplanet with relatively low res- 



- 3 - 



olution, and found that the spin rotation vector of the protoplanet becomes parallel to the 



thors (IMiyoshi et al. I Il999l : iLubow et al. 



1999: 


Kiev et al. 


2001 





2002 



20031 ). However, since the main purpose of these studies was to clarify the planet migration 
process in a large scale (i.e., outside the Hill radius), they did not investigate the region 
near the protoplanet (i.e., inside the Hill radius) with sufficient resolution. Thus, in their 
simulations, we cannot study the gas stream inside the Hill radius. 

To investigate the accretion flow onto the protoplanet system, we need to cover a large 
spatial scale from the region far from the Hill sphere to that in the proximity to the proto- 
planet. For Jupiter, since the Hill radius is th = 744 rp, where Tp is Jovian radius, we have 
to resolve at least ~1000 times different scal es. To cover a large dynamica l range in scales, 
a few authors used the nested-grid method. iD'Angelo et al. I (120021 . l2003l ) investigated the 
relation between the spiral patterns within the Hill radius and migration rate using three- 
dimensional nested-grid code. Although they resolved the region inside the Hill radius, they 
did not investigate the structure in the proximity to the protoplanet because they adopted 
a sink cell at O.lrn that corresponds to ~ 70 rp at Jovian orbit (5 .2 AU). Thus, we can- 
not obs erve a circumplanetary disk at r ^ rn in their calculation. iTanigawa fc Watanabe 
(l2002al ) also investigated the gas flow around a protoplanet using two-dimensional nested-grid 
code. They resolved the region from 12 rn to 0.005 rn. They found that a circumplari etary 
disk with 100 Mm is fo rmed around the protoplanet. ITanigawa fc Watanabd ( j2002bl ) and 
Machida et al. I (l2006d ) investigated the evolution of the protoplanet system using three- 
dimensional nested-grid code. They found that the gas flow pattern in three dimensions is 
qualitatively different from that in two dimensions: the gas is flowing into the protoplanet 
system only in the vertical direction in three-dimensional simulations. 

In the present study, we calculated the evolution of the protoplanet system using three- 
dimensional nested-grid code. We found that after the flow around the protoplanet reaches a 
steady state, the angular momentum accreting onto the protoplanet system is well converged 
regardless of both the cell width and the size of the sink cell region, while the mass accretion 
rate is not well converged. Although we calculated the evolution of the protoplanet system 
with spatial resolution much higher than previous studies, we still need further higher spatial 
resolution to determine the mass accretion rate onto the protoplanet. Thus, in this paper, we 
focus on the gas flow onto and around the protoplanet system and the accretion process of the 
angular momentum, and do not deal with the mass accretion rate (we plan to investigate the 
mass accretion rate with higher spatial resolution using a higher-performance computer in a 
subsequent paper). Note that the speciflc angular momentum accreting onto the protoplanet 
system with fixed protoplanet mass does not strongly depend on the mass accretion rate as 
described in the following sections. 
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The structure of the paper is as follows. The frameworks of our models are given in §2 
and the numerical method is described in §3. The numerical results are presented in §4. §5 
is devoted for discussions. We summarize our conclusions in 56. 



2. MODEL 
2.1. Master Equations 

We consider a local regio n around a protoplanet using shearing sheet model (e.g. 



Goldreich fc Lynden-Belllll965l ). We assume that the temperature is constant and the self- 
gravity of the disk is negligible. The orbit of the protoplanet is assumed to be circular on 
the equatorial plane of the circumstellar disk. 

We set up local rotating Cartesian coordinates with the origin at the protoplanet and 
the X-, y-, and 2;-axis are radial, azimuthal, and vertical direction of the disk, respectively. 
We solve the equations of hydrodynamics without self-gravity: 

|^ + V-(pt;) = 0, (1) 

dv 1 

+ (v\/)v = —VP - V$efr - 2^^p X V, (2) 
ot p 

where p, v, P, ^^s, and ilp are the gas density, velocity, gas pressure, effective potential, 
and Keplerian angular velocity of the protoplanet, respectively. In the above equations, the 
curvature terms are neglected. We adopt an isothermal equation of state, 

P = cIp, (3) 
where Cg is the sound speed. The Keplerian angular velocity of the protoplanet is given by 

where G, M^., and Cp are the gravitational constant, mass of the central star, and orbital 
radius of the protoplanet, respectively. The effective potential $eff is given by 

$eff = -^(3x^-.^)-^, (5) 
where Mp and r are the mass of the p rotoplanet, and the distance from the center of the 



protoplanet (e.g., iMiyoshi et al. Ill999l ). The first term is composed of the gravitational 
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potential of the central star and the centrifugal potential, and higher orders in x, y and z 
are neglected. The second term is the gravitational potential of the protoplanet. Using the 
Hill radius ^ 



equation ([5]) can be rewritten as 



= (7) 



2.2. Circumstellar Disk Model 



Our initial settings are similar to iMiyoshi et al. I (1l999l ) and iMachida et al. I (l2006d ). 

The gas flow has a constant shear in the x-direction as 

«o = (0, -^n^x, 0). (8) 
For hydrostatic equilibrium, the density is given by 



Po 



exp 



(9) 



where uo (= J^^pdz) is the surface density of the unperturbed disk. The scale height h is 
related to the sound speed h = cJQp. 



In the standard solar nebular model (lHayashilll98ll : iHayashi et al. lllQSSi ). the temper- 
ature T, sound speed Cg, and gas density po are given by 

1/4 



280 



lAU/ 



where L and Lq are the protostellar and solar luminosities 

T 



/ kT 



1.9 X 10" 



1/2 /2.34N1/2 



cms 



(10) 



(11) 



10 Ky \ p 

where p = 2.34 is the mean molecular weight of the gas composed mainly of H2 and He, and 



Po = 1.4 X 10 
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lAU/ 



g cm 



(12) 



respectively. When Mc = IMq and L = ILq are adopted, using equations (jlj), ffTOj) . and 
f fTT]) . the scale height h can be described as 

, 5/4 

cm. (13) 



/i = 5.0x10^1 



VlAU 
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2.3. Scaling 

Our basic equations can be normalized by unit time, and unit length, h. The 
density is also scalable in equations (1) and (2) since we neglect the self-gravity of the disk 
in equation (2). We normalize the density by ao/h. Hereafter the normalized quantities are 
expressed with tilde on top, e.g., x = x/h, p = p/{ao/h), t = tflp, etc. The non-dimensional 
unperturbed velocity and density are given by 

v = {0,-h,0), (14) 
Po = ^=exp . (15) 



/2^ V 2, 

Thus, non-dimensional equations corresponding to equations (P), ([2]), ([3]), and ([7]) are 

^ + V-(pi;) = 0, (16) 

Qjj _ 1 _ _ _ _ 

+ {v ■V)i = --VP -V<^es -2z X i, (17) 
ot p 

P = p, (18) 

= -ksx' - ~z') - ^, (19) 
2 r 

where z is a unit vector directed to the 2;-axis. Thus, the gas flow is characterized by 
only one parameter, the non-dimensional Hill radius fn = r-n/h. In this paper, we adopt 
fn = 0.05 — 4.21 (see. Table 1). As functions of the orbital radius and the mass of the central 
star, the parameter fn are related to the actual mass of protoplanet in the unit of Jovian 
mass Mj as 

For example, in the model with th = 1.0, = 5.2 AU and = 1 Mq, the protoplanet 
mass is Mp = 0.4Mj (model M04 in Tabled]). Hereafter, we call model M04 'the fiducial 
model.' For each model, the protoplanet mass for Op = 5.2 AU and Mc = 1 Mq is presented 
in Tabled! In our parameter range, at Jovian orbit {a^ = 5.2 AU), protoplanets have masses 
of 0.05 — 30Mj. We will show our results assuming Op = 5.2 AU and = 1 Mq in the 
following. We will discuss the dependence on the orbital radius Op in §5.21 



In subsequent sections, we use non-dimensional quantities (e.g., p, x, y, and z) when 
we show the structure of the protoplanet system (Figs. 1, 4, 5, 7, 8, 12, and 13). On the 
other hand, to compare physical quantities derived from numerical results with those of 
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the present Jupiter, we use dimensional quantities at Jovian orbit when we show the time 
evolution or radial distribution (Figs. 2, 3, 6, 9, 10, and 11) of the mass M and (specific) 
angular momentum J (j) of the protoplanet system, in which the conversion coefficient from 
physical quantities at Jovian orbit (5.2 AU) into those at any orbit (ap) is described in each 
axis. In addition, for convenience, we redefine the time unit as tp = t/(27r) = l/(2'7rf2p) that 
corresponds to the orbital period of the protoplanet. When ap = 5.2 AU is assumed, tp = 1 
corresponds to 11.86 yr. 



3. NUMERICAL METHOD 

3.1. Nested-Grid Method 

To estimate the angular momentum acquired by a protoplanet system from the proto- 
planetary disk, we need to cover a large dynami c range of spatial scale. Using the nested-grid 



method (for details, see lMachida et al. Il2005al . l2006al ). we cover the region near the proto 



planet by the grids with high spatial resolution, and region remote from the protoplanet by 
the grids with coarse spatial resolution. Each level of rectangular grid has the same number 
of cells (= 32 X 128 x 8), but cell width As{l) depends on the grid level /. The cell width is 
reduced 1/2 with increasing the grid level {I ^ l + l). In a fiducial model, we use 8 grid levels 
(^max = 8). The box size of the coarsest grid / = 1 is chosen to {L^, Ly, L^) = (30, 120, 7.5), 
and that of the finest grid / = 8 is {L^,Ly,L,) = (0.234,0.938,0.059). The cell width in 
the coarsest grid Z = 1 is As = 0.9375, and it decreases with AS = 0.9375/2'^^ as the grid 
level I increases. Thus, the finest grid has AS(8) ~ 7 x 10~^. We assume the fixed boundary 
condition in the x- and i-direction and the periodic boundary condition in the ^-direction. 

In real units, using the standard solar nebular model, the scale height at Jovian orbit 
(op =5.2 AU) is h = 0.27 AU. The computational domain in azimuthal direction corresponds 
to 120 X 0.27 AU =32.4 AU, which is equivalent to the circumference of Jovian orbit around 
the Sun 27rap (32.7 AU). Although we imposed a periodic boundary condition in azimuthal 
direction, this domain size is valid except for ignoring the curvature. Note that the compu- 
tational domain is not necessarily the same as a real circumference of a planet as long as it 
is sufficiently large. To verify our results, in some models, we calculate the evolution of the 
protoplanet system adopting different levels of the finest grid (or different maximum grid 
levels), /max = 5, 6, 7, 9, and 10. In these models, the cell width of /max = 5 is As = 0.23 
(9.375 X 10^'^ cm~^ at Op = 5.2 AU) that corresponds to 33 times Jovian radius, while that 
of /max = 10 is As = 1.83 X IQ-^ (7.32 x 10^ cm^^ at Op = 5.2 AU) that corresponds to 1.02 
time Jovian radius. The maximum grid level /max, and the cell width of the finest grid for 
each model are summarized in Table [TJ 
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3.2. Test Simulation 

At first, we show the evolution of the protoplanet system calculated with Zmax = 8. The 
cell width of the maximum grid level {I = 8) that covers the region in the proximity of the 
protoplanet is As = 7.3 x lO'^ (As = 2.92 x 10^°cm at Op = 5.2 AU). Figure □ shows the 
evolution for model MN04 (fiducial model), in which the protoplanet with 0.4Mj is adopted 
for Op = 5.2 AU. The upper panels in Figure [1] ([T]a - e) show the time sequence of the region 
far from the protoplanet (f > 15) with / = 1, 2, and 3 grid levels, in which three different 
grid levels are superimposed, while the lower panels (Fig. [l]f-j) show the region near the 
protoplanet (f < 3.5) with / = 3, 4, and 5. Each lower panel is 4 times magnification of each 
upper panel. In these panels, the protoplanet is located at the origin (x, y, z) = (0, 0, 0). 
The elapsed time tp is denoted in each upper panel. The central density pc is also denoted 
in each upper panel. 

Figures [T]fl-c and show that the density is enhanced in the narrow band with the 
spiral pattern that is distributed from the upper-left to the lower-right region. The density 
gaps that appears on the right (left) side of the spiral pattern in the region of ^ > {y < 0) 
are also seen in these panels. In Figure [T]c, the density of the spiral pattern around the 
protoplanet is p ~ 2, while that of the gap is p ~ 0.2. Thus, there are a density contrast of 
~ 10 between the spiral patterns and gaps. In addition, the central density increases up to 
p ~ 10^ for tp > 1. Figures HFi-j show a round shape near the protoplanet. Since the Hill 
radius is fn = 1 in this model, the gravity of the protoplanet is predominant in the region 
of f <^ 1. Thus, the central region of f ^ 1 has the round structure. Figure [1] shows that 
the density distributions in panels c-e (or h-j) are similar. Except for the calculation with 
^max = 10, in all models, the structure around the protoplanet har dly changes for tp >, 1 



which seems that the steady state is already achieved for tp > 1. iTanigawa fc Watanabe 



(l2002al ) investigated the evolution of the protoplanet system as the same condition as ours 



but in two dimensions, and showed that the gas stream around the protoplanet is in a steady 



state after a short timescale of tp ~ 1. 



3.3. Convergence Test 

In §3.2[ we showed the evolution of the protoplanet system with the maximum grid 
level /max = 8. In this subsection, to check the convergence of our calculation, we compare 
the evolutions of the protoplanet system with different maximum grid levels (or different cell 
width of the finest grid). Since our purpose is to investigate the angular momentum acquired 
by the protoplanet system, we use the average specific angular momentum as a measure of 
the convergence. As a function of the distance from the protoplanet f, we define the average 
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specific angular momentum jr as 



where the mass 



and angular momentum 



Jr = (21) 



Mr 



Mr 



! Anf'^pdf, (22) 
Jo 



Jr = 47rr p zuV(f, dr, (23) 
Jo 

are integrated from the center f = to distance r. Here, we adopt f = 0.5 fn (i-e., jo.bfu)- 
The dependence of jr on f is discussed in §4.51 We often show the mass and (specific) 
angular momentum of the protoplanet system in real units at = 5.2 AU to compare 
numerical results with present values of gaseous planets. The dimensional values Mr{5.2 AU), 
Jr(5.2 AU), and jr(5.2 AU) at Op = 5.2 AU can be converted into Mr{ap), Jr{ap), and jr{cip) 
at any orbit as 

M,(ap) = M,(5.2AU)(^||^), (24) 
J,(ap) = J,(5.2AU)(^^) , (25) 

and 

jV(ap)=>(5.2AU)(^||^). (26) 

In the following, we show the values at 5.2 AU. When we refer to dimensional physical 
quantities without any mention of orbit, they are the values at Op = 5.2 AU. 

The evolution of jo.srH fo^' models M04L5-M04L9 (see, Table[l]) are shown in FigureO In 
these models, the protoplanet mass is fixed, and only the maximum grid level (or cell width 
of the finest grid) is changed. Figure [2] shows that the average specific angular momentum 
rapidly increases initially, then it saturates at a certain value in each model. Although 
we adopted the same mass of the protoplanet, the saturation levels of the average specific 
angular momenta are different. The average specific angular momentum saturates at jo.srH — 
2 X lO^^cm^ s~^ for model M04L5, while it is saturated at jo.5rH — 7 — 8 x lO^^cm^ s~^ for 
models M04L7, M04, M04L9, and M04L10. Thus, there are a little differences for models 
with /max > 7. For example, model M04L7 has jo.Sm = 6.9 x 10^*^ cm^ s"\ while model M04L9 
has jo.srH = 8-2 x lO^^cm^ s^^ at ip = 10. Therefore, the average specific angular momenta 
are sufficiently converged in /max > 7 or As < 1.4 x 10^^ within a relative error of about 15%. 
In the following, we safely calculate the evolution of the protoplanet system with /max = 8 
of the maximum grid level. 
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3.4. Models with and without Sink Cells 

When we adopt the maximum grid level /max = 8, the cell width is AS = 7.3 x 10^^. 
In real units, when the protoplanet is located at 5.2 AU, the cell width corresponds to 
As = 2.9 X 10^°cm whose size is 4.1 times the Jovian radius. Interior to the gas gian t 



planet, there exists a solid core with ~ 10 in mass and ~ 10^ cm in size ( lMizundll980l ). 
The size of the solid core is much smaller than the cell width adopted in our calculation. 
To investigate the evolution of a protoplanet through gas accretion, in principle, we need 
to resolve a solid core with sufficiently small cell size (As ^ 10^ cm). However, since 
our purpose is to investigate the angular momentum flowing into the protoplanet system, 
we do not always n eed t o resolve a central solid core and a protoplanet. As discussed 



m 



D'Angelo et al. I (120021 ) . the gas around the central region makes an artificial pressure 
gradient force, which may affect the gas accretion onto the protoplanet system. To check 
this, we also calculated the evolution of the protoplanet system adopting the sink cell in 
some models. We parameterized the size of the sink: fsink = 0.01 and 0.03 (models M04S01 
and M04S03). These models are also summarized in Tabled! 

In real units, the radius of the sink in model M04S01 is 4.0 x 10^'^ cm (~ 5.6 Jovian 
radius), while that in model M04S03 is 1.2 x 10^^ cm ( ~ 17 Jovian radius). During the 
calculation, we remove the gas from the region inside the sink radius in each time step, 
and integrate the removed mass and angular momentum that are assumed as the mass and 
angular momentum of the protoplanet system. 

After the steady states are achieved at tp ~ 20, we estimate the average specific angular 
momenta jr as a function of f. Figure [3] shows the distribution of the average specific angular 
momentum jj. against the distance from the protoplanet r in real units for each model. The 
solid line represents jr without a sink cell, while the other lines represent those with sink 
cells. Inside the sink radius, the average specific angular momentum is assumed as a constant 
value. In Figure [3], the vertical dotted line indicates the Hill radius. 

In Figure [3] the solid line shows that jr increases from the center to a peak around 
the Hill radius, and drops sharply just outside the Hill radius. The drop indicates that 
the angular momentum becomes negative at r > rn. Thus, the rotational d irection turns 



around between the region in side and outside the Hill radius. A s shown in ISekiya et al 



(119871 ) , iMiyoshi et al. I (119991 ) , and iTanigawa fc Watanabd (l2002al ) , the protoplanet formed 



by the gas accretion in the protoplanetary disk has a prograde spin, and thus it has a positive 
(specific) angular momentum. On the other hand, gas far outside of Hill sphere seems to 
rotate retrogradely against the protoplanet because it rotates with nearly Keplerian velocity 
with respect to the central star [v = —(3/2) QpX, as shown in eq. ([HD]. As a result, gas inside 
the Hill radius has a positive angular momentum, while that outside the Hill radius has a 
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negative angular momentum. 

Figure E] shows that jr in models with the sink cell depends on the radius of the sink 
in the region r <^ th- However, they do not depend sensitively on the size of the sink in 
the region of r > 0.5 th- For example, at r = rn, model without the sink (model M04) has 
jj.^ = 5.3 X lO^'^cm^ s"\ while model M04S03 has j^^ = 6.4 x lO^'^cm^ s"^ Thus, difference of 
the average specific angular momentum among these models at the Hill radius is ~ 12%. This 
difference decreases with the sink radius. In this study, we focus on the angular momentum of 
the protoplanet system, not the planet itself. As shown in Figure [31 the angular momentum 
acquired by the protoplanet system (planet + circumplanetary disk) extends up to the Hill 
radius, in which almost all the angular momentum distributed in the region of r ^ rp or 
r ^ Tsink- In the following, we calculate the evolution of the protoplanet system without the 
sink cell. 



4. RESULTS 

4.1. Typical Gas Flow 

Figure H] shows the gas structure around the Hill sphere for model M04 at tp = 20. 
Figure m left panel shows the structure on the cross section in the z = plane, in which red 
lines indicate the streamlines. In this panel, the gas enters / = 4 grid from upper y boundary 
for X > (from lower y boundary for 5; < 0) and goes downward (upward for x < 0) according 
to the Keplerian shear motion. The shocks (crowded contours near the Hill radius) are seen 
in the upper right and lower left region from the protoplanet. In this model, the shock 
front almost corresponds to the Hill radius (Fig. H] left panel) . When the gas approaches 
the protoplane t , the streamlines are bent by the gravity of the protoplanet. According to 



Miyoshi et al. I (119991 ) . the gas flow is divided into three region: the pass-by region (|x| > rn) 



the horseshoe regi on (x < rn, and 1 7/1 > rn), and the planet atmosphere region (f < rn) 



Note that although iMiyoshi et al. I (119991 ) classified the flow pattern in their two dimensional 
calculation, their classification is useful for global flow pattern in three dimensions. In the 
pass-by region, the flow is first attracted toward the protoplanet, and then causes a shock 
after passing by the protoplanet. At the shock front, the density reaches a local peak and the 
streamlines bend suddenly. On the other hand, the gas entering the horseshoe region turns 
round by the Coriolis force and goes back. The outermost streamlines in the horseshoe region 
(i.e., the streamlines passing very close to the protoplanet) pass through the shock front, 
while the gas on the streamlines far from the protoplanet does not experience the shock. In 
the atmospheric region, the gas is bound by the protoplanet and forms a circumplanetary disk 
that revolves circularly around the protoplanet in the prograde (counterclockwise) direction. 
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Although the streamhnes on z = plane (Fig. H] left panel) are similar to those in recent 



two dimensional calculations (e.g.. lLubow et al. Ill999l : iTanigawa &: Watanabdl2002a] ). there 



are important differences. In two-dimensio nal calculation s, a pa rt of the gas near the Hill 
sphere can accrete onto the protoplanet. iLubow et al. I (Il999l ) showed that gas only in 
a narrow band that distributed from the lower left to the upper right region against the 
protoplanet for y < spirals inward toward the protoplanet passing throug h the shocked 



regio n and finally accretes onto the protoplanet (for details, see Figs. 4, and 8 of lLubow et al. 



19991 ). On the other hand, in our three-dimensional calculation, gas only flows out from the 
Hill sphere and thus does not accrete onto the protoplanet on the midplane. Figure H] left 
panel shows that although gas flows into the Hill sphere, a part of the gas flows out from 
the central region. Figure H] right panel is three-dimensional view at the same epoch as the 
left panel. In this panel, only the streamlines flowing into the high-density region of f <^ fn 
are drawn for z > which are inversely integrated from the high-density region. This panel 
clearly shows the gas flowing into the protoplanet in the vertical direction. 

To investigate gas flowing into the protoplanet system in detail, in Figure [5l we plot 
the streamlines at the same epoch as Figure H] with different grid levels (/=3, 5, and 7). In 
this figure, each upper panel shows three-dimensional view, while each lower panel shows the 
structure on the cross section in the y = plane. Note that, in lower panels, the streamlines 
are projected onto the ^ = plane. Figure [5]a shows only the streamlines in a narrow bundle 
flowing into t he protoplanet syste r n. This feature is similar to tha t shown in two-dimensional 
calculations (ILubow et al. Ill999l : iTanigawa fc Watanabdl2002al ). However, the streamlines 
in Figure Eh indicate that gas rises upward near the shock front and then falls into the central 
region in the vertical direction. Gas flowing in the vertical direction spirals into the inner 
region (Fig.[5]c). In this process, vortices appear as shown in Figure [5](i. As shown in Figure H] 
left panel, also in Figure [5]d, gas is flowing out from the central region on 5 = plane. Gas 
in the proximity of the protoplanet rotates circularly in the prograde direction as shown in 
Figure Oe. Figure [S]f shows that a part of the gas flowing into the upper boundary of / = 7 
grid level contributes to the disk formation around the protoplanet, while a remainder is 
bent and flows out from the central region. 

When we look down the protoplanetary disk from the above along the 5-axis, streamlines 
may seem to be almost the same as those in two-dimensional calculations. However, gas 
moves also in the vertical direction: streamlines go upward near the shock front (f ~ fn), 
and vertically falls into the central region at f ^ fn. This feature of streamlines is also 
seen in models M04S01 and M04S03, in which the sink cell is adopted. Thus, different 
features of streamlines in two- and three-dimensional calculations are not caused by the 
effect of the pressure gradient force, but caused by the dimensions (because the same feature 
appears in both models with and without the sink). This flow pattern is also seen in other 
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three-dimensional calculation (iTanigawa &: Watanabell2002bl : lD'Angelo et al. II2003I ). In two- 
dimensional calculation, since the vertical motion is restricted, the flow pattern is different 
from that in three din aensions. This difference affects the accretion rate onto the protoplanet 
and migration rate. iD'Angelo et al. I (120031 ) showed that the migration rate is different 
between two- and three-dimensional calculations. In the present study, however, since we 
focus on the angular momentum of a protoplanet system, we do not discuss them any more. 
We will discuss the accretion and migration rate in the subsequent papers. 

Finally, we comment on the circumplanetary disk. In Figure [5] lower panels, the green 
surface (i.e., iso-density surface) indicates the high-density structure around the protoplanet. 
These panels show the disk-like structure in the proximity of the protoplanet, and the disk 
becomes thinner as it approaches the protoplanet (i.e., the origin). We discuss the circum- 
planetary disk in §5.31 



4.2. Dependence on Protoplanet Mass 



We have shown the evolution of the protoplanet system for model with 0.4Mj in §4.11 
In this subsection, we investigate the evolution of the protoplanet system with different 
protoplanet masses. Figure [6] upper panels show the accumulated masses (eq. [22]) within 
f < 0.1 (Mq.i, Fig. [H]a) and f < 0.05 (M0.05, Fig. ^) against the elapsed time for different 
models, while Figure [6] middle panels show the corresponding angular momenta (eq. [23l) in 
the same regions (Fig.[6]cfor Jq.i, and Fig.[6]dfor Jq.os)- In Figure [6l both masses and angular 
momenta for all models increase rapidly for tp < 0.1. This is because the protoplanet with 
mass of 0.05 — 0.6Mj suddenly appears in the protoplanetary disk at tp = 0. However, 
this rapid growth phase (tp < 0.1) is not real, because the gas is considered to begin to 
accrete onto the protoplanet wheri the mass of the so lid core exceeds ~ 10 (IMizuno 
1980l : iBodenheimer fc Pollacklll986l : llkoma et al. Il2000l ). The growth rates of the mass and 



angular momentum begin to decrease at tp ~ 0.1 in all models, then both masses and 
angular momen ta increase with an aln i ost con stant rate until the end of the calculation 
(0.1 < tp < 20). ITanigawa fc Watanabd (l2002al ) calculated the mass accretion rate onto the 
protoplanet as 



M = 8.0 X lO^^Mff 



lOMff 



1.3 



(27) 



Thus, the growth time Xgrow = M/M is Xgrow = 1000 — 430 yr (i.e., tp = 88 — 36). In our calcu- 
lation, we continue to calculate the evolution of the protoplanet system for t ~ 230 yr (tp ~ 
20) by fixing the planet mass. We think that this treatment is not problematic, since the 
growth timescale is longer than our calculation time. Note that since ITanigawa &: Watanabe 
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( l2002al ) calculated the evolution of the protoplanet system in two dimensions, the growth 



rate might be different from that in three-dimensional calculations. 

Figure M lower panels show the evolution of the average specific angular momentum 
in the region of f < 0.1 (Fig. [6]e for jo.i) and f < 0.05 (Fig. [6]f for jo.os)- The masses 
and angular momenta flowing into the protoplanet system increase with a constant rate for 
tp < 0.1, while the average specific angular momenta are saturated at certain values for 
ip > 0.1. This saturation means that flow around the protoplanet is in the steady state. 
Figures [6]e and / also indicate that the average specific angular momentum brought into the 
protoplanet system increases with the mass of the protoplanet. In Figures [6]e and /, the 
average specific angular momenta in the region of f < 0.1 is larger than those in the region 
of f < 0.05 indicating that the protoplanet system has a larger average specific angular 
momentum in the more distant place from the protoplanet. We will investigate the angular 
momentum acquired in the protoplanet system in §4.51 



4.3. Gas Structure around a Protoplanet 



Figures [7] and [8] show the density distributions (upper panels) and Jacobi energy con- 
tours (lower panels) on the cross section in the z = (Fig. [7]) and y = plane (Fig. [8]) 
around the Hill sphere after the steady state is achieved (tp ~ 20) for models M02 (left 
panels), M04 (middle panels), and M06 (right panels). The white-dotted line in each panel 
indicates the Hill radius fn. In each upper panel, the shock appears from the upper left 
to lower right near the Hill radius. These shocks are frequently seen in similar calculations 



(ISekiva et al. 



1987 



Mivoshi et al. Ill999l : iLubow et al. Ill999l : iTanieawa fc Watanabell2002a 



D'Angelo et al. II2002I ). In Figure [7] upper panels, the round structures are seen in the prox- 
imity to the protoplanet (f -C th), while the ellipsoidal structures are seen in the regions 
of f ~ th. This is because gas distributed near the protoplanet is more strongly bound by 
the protoplanet. Figure [8] upper panels show that contours of the central region sags in the 
center of a concave structure, and thin disks are formed around the protoplanet (f ^ fn). 
In addition, the butterfly-like structure is also seen inside the Hill radius in Figure [8] upper 
panels. These structures are considered to be formed by the rapid rotation of the central 
circumplanetary disk: s imilar structure is seen in a rapidly rotating protostar (e.g., Fig.l of 
Saigo fc Tomisakall2006l ). We will discuss the disk structure in §5.31 



In contrast to celestial mechanics, it is difficult to find fluid elements bound by the 
protoplanet because thermal energy is important in addition to the gravitational and kinetic 
energies. To discern gas bound by the gravity of the protopla net, we use the Jacobi energy as 
an indicator (ICanup fc Esposito Ill995l : iKokubo et al. Il2000l ). In our unit, the Jacobi energy 
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is given by 

1 /^~2 -^2 ■^2\ 3^2 1-2 3f?r 
^^ = 2r +^ +'')~2 +2 ^ ^ 

In equation fl28l) . the first term is the kinetic energy, the second and third terms are the tidal 
energy of the central star, and the forth term is the gravitational energy of the protoplanet. 
The Jacobi energy is the conserved quantity in a rotating system. In the case of fluid, the 
Jacobi energy is not strictly appropriate because the thermal energy is ignored. However, 
we can use this for rough estimation of gas bound by the protoplanet. We determine fluid 
elements bound by the protoplanet from the contour of the Jacobi energy in the lower panels 
in Figures [7] and [HI 

Fluid elements with lower Jacobi energy are strongly bound by the protoplanet as de- 
scribed in equation (1251) . Lower panels of Figures [7] and M show that inside the region of 
f < 0.5 fn, each contour has a closed ellipse. Since fluid elements move on this closed orbit, 
it is considered that gas distributed in the region of f < 0.5fH is bound by the protoplanet 
when the thermal effect can be ignored. We discuss the thermal effect in §4.51 On the other 
hand, although fluid elements exist inside the Hill radius, outside 0.5 fn ^ ^ f^n, they 
are not bound by the protoplanet, because the contours of the Jacobi energy straddle the 
Hill radius. For example, in lower left panel of Figure [TJ the contour of Ej = —2.5 straddle 
the Hill radius, and thus fluid elements with this Jacobi energy freely move on this contour. 
Thus, although these elements transiently stay inside the Hill radius, they flow out from the 
Hill sphere. In summary. Figures [7| and [S] suggest that fluid elements inside f < 0.5 fn are 
bound by the protoplanet. 

4.4. Angular Momentum of a Protoplanet System 

Figure [9] upper panel shows the accumulated mass Mr for different models after the 
steady state is achieved (tp ~ 20). In this panel, a thin solid line indicates the initial value 
(or the value of the protoplanetary disk), and the circle is the Hill radius rn for each model. 
The accumulated mass in any model is larger than the initial value, because gas flows into 
the Hill sphere. This panel indicates that the massive protoplanet has a massive envelope. 
Since these mass distributions are in a steady state at the fixed mass of the protoplanet, it 
can be considered that different curves correspond to snapshots at different evolution phases. 
Namely, a gas envelope increases its mass with time (or the protoplanet mass). Outside the 
Hill radius, the accumulated mass in each model converges to the initial value indicating 
that the mass distribution for r ^ rn does not change from the initial state because the 
influence of the protoplanet is small. 
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Figure [9] lower panel shows the absolute value of the angular momentum | | . There 
are spikes in all models, at which the sign of the angular momentum is reversed. As shown 
in §3.41 the angular momentum has a positive sign around the protoplanet (r <^ th), while 
it becomes negative outside the Hill radius (r ^ th). The sign of the angular momentum 
is reversed outside the Hill radius for models MOl, M02, M04, and M06, while it is reversed 
inside the Hill radius for model MOOS. This is because the protoplanet system for model MOOS 
does not acquire a sufficient mass and angular momentum owing to the shallow gravitational 
potential and relatively large thermal pressure (for details, see §4. SI) . Figure [H] upper panel 
shows that the protoplanet system has the envelope mass of only M ~ 0.06 for model 
MOOS, in which the protoplanet mass is SMj. Except for model MOOS, the angular momenta 
gradually decrease after they reach their peak around r ~ th, then it becomes negative at 
r ~ 1 — 2rH. Thus, the angular momenta bound by the protoplanet system are limited in the 
region of r < 2rH at the maximum. Figure lower panel shows that the angular momentum 
keeps an almost constant value around the Hill radius (O.Srn r < It'h), until they are 
reversed. This means that the angular momentum with plus sign and that with minus sign 
are mixed in this region, as shown in Figure H] left panel. As a result, wherever we estimate 
the angular momentum in the range of O.Srn t < I^h, we can obtain almost the same 
values of the angular momentum. 

The distributions of the average specific angular momentum jV- for different models are 
shown in Figure [TOl The circles in this figure mean the Hill radii rn. The crosses indicate the 
Jacobi radii rj inside which gas is considered to be bound by the gravity of the protoplanet. 
We determine the Jacobi radii from the contour of the Jacobi energies as in Figures [7| and [S] 
lower panels. The Jacobi radii in all models are distributed in the range of rj ^ O.S — Irn. 
Figure [10] indicates that a more massive protoplanet has an envelope with larger amount of 
the specific angular momentum. Thus, the specific angular momentum accreting onto the 
protoplanet system increases as the protoplanet mass increases. The rapid drops at large 
radii indicate the reverse of the rotation axis as shown in Figures [3] and O In Figure (TUl in 
model MOl (M = O.lMj), > at r = rn is twice of that at r = rj, while, in models M02, M04, 
M06, Ml, and M3 (M > 0.2Mj), there are little differences between the average specific 
angular momentum derived from the Hill radius jVjj and that derived from Jacobi radius jV-i- 
Thus, we can safely estimate the average specific angular momentum using either the Hill 
radius or Jacobi radius for models with M > 0.2Mj. 
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4.5. Evolution of the Specific Angular Momentum 



To properly calculate the angular momentum of the protoplanet system, we have to 
calculate the planetary growth from the solid core with ~ 10 to the present mass. How- 
ever, it takes huge computation time to calculate all evolution phases. Thus, we estimate 
the angular momentum of the protoplanet system according to the following procedure: (i) 
we calculate the average specific angular momenta of the gas fiowing into the protoplanet 
system at the fixed masses of the protoplanet (i.e., under the same parameter fn) in mod- 
els changing the mass of the protoplanet, then (ii) derive the relation between the average 
specific angular momentum and mass of the protoplanet, and describe it as a function of 
the protoplanet mass, and (iii) estimate the angular momentum of the protoplanet system 
integrating the average specific angular momentum by mass up to the present value of gas 
giant planets. 

At first, we analytically estimate the specific angular momentum of the protoplanet 
system, then compare it with numerical results. We assume that the gas that overcomes the 
Hill potential fiows into the Hill sphere (or a protoplanet system) with the Kepler velocity of 
the protoplanet at the Hill radius rn. When the protoplanet mass is Mp, the Kepler velocity 
at th is given by 



Note that is constant when Op is fixed. The specific angular momentum can be written 

as 



Thus, the specific angular momentum is proportional to oc r^. In addition, when the mass 
of the central star is fixed, the specific angular momentum is proportional to 



Thus, the specific angular momentum increases with 2/3 power of the protoplanet mass. 

Figure [TT] shows the average specific angular momentum derived from all models against 
the protoplanet mass. The average specific angular momenta are estimated in the region 
of r < rn (+) , r < O.Srn (A), r < 0.1 rn (□), and r < rj (o). Although we fixed the 
protoplanet mass in each model, we can consider the horizontal axis in Figure [11] as time 
sequence of the protoplanet system. Figure [TT] clearly shows that more massive protoplanet 
can acquire the envelope with larger average specific angular momentum. Although the 
average specific angular momenta for r < O.Srn differ from those for r < rn for models with 
Mp < 0.2Mj, there are little difference for models with Mp > 0.2Mj. In addition, when we 
adopt the average specific angular momenta in the region of r < 0.1 rn, we underestimate 




(29) 



(30) 



J oc M2/3. 

J p 



(31) 
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them for any model as shown in Figure [TTl Therefore, we can properly estimate the average 
specific angular momentum of the protoplanet system in any region of r < 0.5 — Ith for 
models with Mp > 0.2Mj. 

In Figure [TTl the red and blue lines are fitting formulae of the evolution of the average 
specific angular momentum as a function of the protoplanet mass for Mp < Mj (low mass, 
jiin; red line), and Mp > Mj (high mass, jhm: blue hne). They are given by 



Figure [TT] shows that the evolution of the average specific angular momentum for Mp > IMj 
is well described by j oc Mp'''^ which corresponds to the equation ( l3Ti) . On the other hand, 
the growth rate of the average specific angular momentum for Mp < IMj is larger than 
j oc Mp''^, and can be fitted by j oc Mp. We ignored the thermal effect when we derive 
equation (l3T|l . When Mp is small, the gas flowing into the protoplanet system is affected 
relatively strongly by the thermal pressure. We can estimate the mass at which the gravity 
dominates the thermal pressure force from the balance between the thermal pressure gradient 
and gravitational forces. Near the Hill radius, the thermal pressure gradient force is more 
dominant than the gravity when the Kepler speed is slower than the sound speed (i.e., 
< Cs). On the other hand, when vk > Cs, the gas flow is controlled mainly by the gravity 
of the protoplanet even near the Hill radius. Using equation fl2^ . fx = Cg is realized when 
'^H = Cs/{y/^^p) which corresponds to M = 0.08Mj of the protoplanet mass at Op = 5.2 AU. 
Thus, the gas flow is largely affected by the thermal pressure for Mp <C 0.08Mj, while it is 
not so affected by the thermal pressure for Mp ^ 0.08Mj. In Figure [TT| the evolution of the 
angular momentum for Mp > IMj well corresponds to the analytical solution {j oc Mp*^^) 
indicating that the thermal effect is negligible for Mp > IMj. However, for Mp < IMj, 
the average specific angular momentum is smaller than jhm (blue line). Thus, when the 
protoplanet is Mp < IMj, the thermal effect is not negligible for the acquisition process of 
the angular momentum because the gravitational potential is relatively shallow. Thus, the 
thermal pressure seems to remain important for M < IMj, although vk = Cg realized at 
M = 0.08Mj. 

To verify the relation of the thermal and gravitational effects, in Figure [121 we plot the 
ratio of the azimuthal to Kepler velocity (f^/fx) around the protoplanet on the midplane 
for models M02, M04, and M06. In this figure, closed contours inside the Hill radius indicate 
that gas revolves around the protoplanet. For example, the closed contour of v^/vk = 0.5 in 




for Mp < IM 



(32) 



and 




(33) 
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Figure [T2l means that gas rotates along the contour with 50% of the Kepler velocity. The black 
circles represent contours of = Cg, inside which gas revolves around the protoplanet with 
supersonic velocity. Thus, the gravitational force of the protoplanet is dominant inside the 
black circles, while the thermal pressure gradient force is dominant outside the black circles. 
Figure [12] shows that the radius of the f ^ = Cg contour increases with the protoplanet mass 
indicating that the region that dominated by gravity of the protoplanet extends with the 
protoplanet mass. As the protoplanet mass increases, the flow speed inside the Hill radius 
approaches to the Kepler velocity, and the region dominated by the gravity of the protoplanet 
spreads outward. In this way, since the thermal effect decreases as the protoplanet mass 
increases, the growth rate of the average specific angular momentum approaches to the 
analytical solution (eg. [311). 



5. DISCUSSION 

5.1. Effect of Thermal Pressure 

In this paper, we investigated the evolution of the protoplanet system under the isother- 
mal approximation, which might not be problematic in the circumplanetary disk, but may 
not be valid in the very proximity to the protoplanet. To investigate the thermal effect 
around the protoplanet, we adopted the sink cell in some models. As shown in §3.4[ the 
cell width in the fiducial model having /max = 8 is about 4 Jovian radius. Thus, gas inside 
r < 4rjup has the same thermal energy. If the protoplanet has almost the same size as the 
present gas giant planet, the thermal energy around the protoplanet may be overestimated 
in model without the sink cell, because an actual protoplanet is embedded in the small part 
of the innermost cell. On the other hand, when we adopt the sink cell, the thermal energy 
around the protoplanet is underestimated, because the thermal energy is artificially removed 
from the sink cell. However, as shown in §3.4[ when the sink radius is much smaller than the 
Hill radius (rsink "C l/lOOrn), there are little differences in the angular momentum acquired 
by the protoplanet system between models with and without the sink cell. This is because 
a large part of the angular momentum of the protoplanet system is distributed around the 
Hill radius (0.5 ^ t < th) as shown in §4.41 The Jovian radius is much smaller than the 
Hill radius (rn = 744rjup). In addition, the angular momentum flowing into the protoplanet 
system is determined by the shearing motion in the region of r ~ th (see 14.51) . Therefore, 
under the assumption that the protoplanet is smaller than the innermost cell width, it is 
expected that the thermal effect from the protoplanet is sufficiently small for the acquisition 
process of the angular momentum when we are using sufficient smaller cells than the Hill 
radius. 
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Mizund (Il980l ) , iBodenheimer k. Pollackl (119861 ) , and llkoma et al. I ( I2OOOI ) sugges ted that 
gaseou s protoplanet has a large envelope with high temperature. For example, iMizuno 
( llQSOl ) showed that, for Jovian case, gas distributed in the range of r > 2 x 10^^ cm behaves 
isothermally, while that distributed in the range of r < 2 x 10^^ cm b ehaves adiabatically 
in their spherically symmetric calculations (see Fig. 4 of iMizund Il980l ). Thus, gas in the 
range of r < 2 x 10^^ cm (~ 30 rjup) has higher temperature than the ambient medium. 
However, it is expected that the gaseous envelope hardly affect the angular momentum of 
the protoplanet system, because the size of the envelope (~ 30 rjup) is much smaller than the 
Hill radius (rn = 744 rjup). A large part of the angular momentum is distributed in the range 
of O.Srn T < th, as shown in §4.4[ On the other hand, the size of the circumplanetary 
disk may be affected by the thermal envelope. Thus, when we investigate the formation 
of the circumplanetary disk, we have to include the realistic thermal evolution around the 
protoplanet. However, it is difficult to study the thermal evolution around the protoplanet, 
because we need to solve the radiation hydrodynamics in three dimensions. In a subsequent 
paper, we will investigate the effect of the thermal envelope under simple assumptions. 



5.2. Angular Momentum for Jupiter and Saturn 

In §4.51 we fixed the physical quantities as those at Op = 5.2 AU. However, in our 
calculation, since we use th e dimensionless quant ities, we can rescale those at any orbits a^. 



Under the standard model (IHayashi et al. Ill985l ). we can generalize equation (132|) as 



Jlr 



7.8 X 10- (MA (-^)'"cmH-' = 2.5 X 10- (^) (-^)"%n,= .-'. (34) 



We assume that the protoplanet system acquires the gas with the average specific angu- 
lar momentum of equation (1341) when the protoplanet has a mass of Mp < IMj. Thus, 
to estimate the angular momentum of the protoplanet system, we need to integrate equa- 
tion fl34|) by the mass until the present values of gas giant planets. For example, the angular 
momentum in our model for Jupiter is 

Jj = / jim(5.2 AU) dM = 1.3 X lO^^m^ s'S (35) 



which is about 30 times larger than that of present Jupiter (4.14 x lO^^gcm^s ^). On the 
other hand, the angular momentum in our model for Saturn is 

/Ms 
jim(9.6AU) dM = 3.6 x 10^^ gcm2s~\ (36) 



- 21 - 



which is 50 times larger than that of present Saturn (7.2 x lO'^^gcm^ s~^), where Ms = 
95.16 is the Saturnian mass. Note that the orbital angular momenta of the present 
Jovian and Saturnian satellites can be ignored in the protoplanet system because they are 
considerably smaller than the Jovian and Saturnian spin angular momenta. 

The above estimation corresponds to the angular momentum of proto planet-disk sys- 
tem, which includes the central planet and protoplanetary disk. We showed that the angular 
momentum of the protoplanet system is 30 — 50 times larger than present spin of the gaseous 
planet. The angular momenta flowing into the Hill sphere are distributed into the spin of the 
protoplanet and orbital motion of the circumplanetary disk. It is expected that a large frac- 
tion of the total angular momentum contributes to the formation of the circumplanetary disk 
and the residual contributes to the spin of the planet. Although the fraction of the angular 
momentum of the circumplanetary disk to the spin of protoplanet is not correctly estimated 
in our calculation, a part of the angular momentum is certainly distributed into the disk. 
Thus, the angular momentum transfer and dissipation mechanism for the circumplanetary 



disk a re necessary to follow further evolution of the protoplanet system. iTakata fc Stevenson 



( 1l996l ) proposed the despin mechanism of the protoplanet by planetary dipole magnetic field, 
in which an initially rapidly rotating protoplanet can be spun down to the present value by 
the magnetic interaction between the protoplanet and circumplanetary disk. 



5.3. Disk formation and Implication for Satellite Formation 

When the circumplanetary disk is formed around the protoplanet, it is possible to 
form satellites in the disk. While there are many scenarios for the satellite formation 



(e.g., IStevenson et al. Ill986l ). regular satellites around gas giant planets are supposed to 
be formed in the gaseous disk as for the plane t formation in the protoplanetary disk (e.g., 
Korycansky et al. Ill99ll : ICanup fc Wardll2002l ). Observations showed that the regular satel- 
lites around Jupiter and Saturn are distributed only in the close vicinity of the planet 
(r < 50 Tp), and are on prograde orbits near the equatorial plane. Thus, it is expected 
that these regular satellites formed in the circumplanetary disk. 

We discussed the angular momentum of the protoplanet system in §4.51 As the angular 
momentum flowing into the Hill sphere is brought into both the protoplanet and circumplan- 
etary disk, we cannot estimate the fraction of angular momentum for the circumplanetary 
disk. When we assume that the centrifugal force is balanced with the gravity of the proto- 
planet, we can derive the centrifugal radius as 



GMr. 



(37) 
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To quantify the disk size, we adopt the specific angular momentum in equation (1371) as those 
in equation flMl) . The centrifugal radius for the proto- Jovian disk with Mp = IMj and 
Op = 5.2 AU is 

rcf,j = 1.5 X 10^^ cm, (38) 

which is 22 times as large as Jovian radius (rjup = 7.1 x lO^cm"^). The Galilean satellites, 
which are regular satellites, are distributed in the range of 6rjup < r < 27rjup, which is 
consistent with the centrifugal radius derived from our calculation. Note that the disk size is 
expected to be larger than equation fl55]) owing to the thermal effect around the protoplanet. 
In the same way, we estimate the centrifugal radius of proto-Saturnian disk as 

rcf,s = 4.1 X 10^^ cm, (39) 

which is 68 times as large as Saturnian radius (rsat = 6.0 x 10^ cm). The Saturnian repre- 
sentative regular satellites (Mimas, Enceladus, Tethys, Dione, Rhea, Titan, Hyperion, and 
lapetus) are distributed in the range of 3rsat ^ ^ ^ SOrsat, which is also corresponding to 
the centrifugal radius of our result. 

Figure US] shows the structure around the protoplanet at the end of the calculation 
(tp ~ 20) for model Ml, in which the protoplanet has mass of IMj. Each right panel is 4 
times magnification of each left panel. Figures [T3la. 6, and c show the shock fronts outside the 
Hill radius, and a thick disk between the shock front and the Hill sphere. In the proximity 
to the protoplanet, a concave structure is seen in Figures [T3]e and /. The contours in these 
figures rapidly drop around the rotation axis (i.e., z-axis). Thus, even inside the Hill sphere, 
except for the central region, the disk has a thick torus-like structure. Figure [T5] q. h, and i 
show the structure in the very proximity to the protoplanet. On the midplane, the elliptical 
structure is seen in a large scale (Fig. [T3](i). while an almost round shape is seen in the 
proximity to the protoplanet (Fig. [T^) because the force acting on gas in this region is 
dominated by the gravity of planet. The contours in Figures [T^ and i show a very thin 
disk in the region of f < 0.1 (or r < 50rjup). Thus, the regular satellites might be formed 
in these circumplanetary disks. However, to study the satellite formation in more detail, we 
need more realistic calculations. 



6. SUMMARY 

To study the gas flow pattern and angular momentum accretion onto the protoplanet 
system, we have calculated the evolution of the protoplanet system in the circumstellar 
disk. Firstly, we have investigated the dependence of the angular momentum accreting onto 
the protoplanet system on spatial resolution with different cell widths, and confirmed its 
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convergence for the cell width much smaller than the Hill radius. Next, adopting the sink 
cell whose size is comparable to or slightly larger than the radius of the present gaseous 
planets in the solar system, we have checked that, in both models with and without the sink, 
the thermal effect around the protoplanet barely affects the gas ffow pattern and angular 
momentum of the protoplanet system. Thirdly, with sufficiently high spatial resolution, we 
have calculated the evolution of the protoplanet system with the protoplanet mass in the 
range of 0.05Mj < Mp < 30Mj, where Mp and Mj are the protoplanet and Jovian masses, 
respectively. The following results are obtained: 

• The gas flow pattern in three dimensions is qualitatively different from that in two 
dimensions: the gas is flowing onto the protoplanet system mainly in the vertical 
direction in three-dimensional simulations. 

• The specific angular momentum increases as j oc Afp when the protoplanet mass is 
Mp < IMj, while it increases as j oc Mp^^ when Mp > IMj. 

• The angular momentum of the protoplanet system is 30 — 50 times larger than the 
present spin of the gaseous planets in the solar system. 

• A thin disk is formed only in the region of r < 20 — 60 Tp, where Vp is the radius of 
the planets. This location agree with the orbital radii of regular sateUites around the 
present gaseous planets in the solar system. 

These conclusions imply that the protoplanet system can acquire sufficient angular mo- 
mentum from the circumstellar disk, and dcspin mechanisms of the protoplanet system arc 
necessary to realize the present spin of the gaseous planets in the solar system. We are 
planning higher-resolution simulations to further investigate the evolution of the angular 
momentum of the gas giant planets in near future. 

We have greatly benefited from discussion with H. Tanaka, and T. Muto. We also thank 
T. Hanawa for contribution to the nested-grid code. Numerical calculations were carried out 
with Fujitsu VPP5000 at the Center for Computational Astrophysics, the National Astro- 
nomical Observatory of Japan. SI is grateful for the hospitahty of KITP and interactions 
with the participants of the program "Star Formation through Cosmic Time" . This work is 
supported by the Grants-in-Aid from MEXT (15740118, 16077202, 16740115, 18740104). 
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Table 1: Model parameters 



Model 




Mp (5.2 AU) 


^max 


As (10"^) 


r-sink (10-2) 


MOOS 


0.5 


0.05 


8 


7.3 




MOl 


0.63 


0.1 


8 


7.3 




M02 


0.8 


0.2 


8 


7.3 




M04 


1.0 


0.4 


8 


7.3 




M06 


1.15 


0.6 


8 


7.3 




Ml 


1.36 


1 


8 


7.3 




M3 


1.95 


3 


8 


7.3 




MIO 


2.91 


10 


8 


7.3 




M30 


4.21 


30 


8 


7.3 




M04L5 


1.0 


0.4 


5 


5.9 




M04L6 


1.0 


0.4 


6 


2.9 




M04L7 


1.0 


0.4 


7 


1.5 




M04L9 


1.0 


0.4 


9 


3.7 




M04L10 


1.0 


0.4 


10 


1.8 




M04S01 


1.0 


0.4 


8 


7.3 


1 


M04S03 


1.0 


0.4 


8 


7.3 


3 



"in unit of Jupiter mass Mj 
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Fig. 1. — Time sequence for model M04. Density {color scale) and velocity distribution 
(arrows) on the cross section in the z = plane are plotted. Each lower panel (/ = 3) is four 
times spatial magnification of each upper panel (/ = !). Three levels of grids are shown in 
each upper (/ =1, 2, 3) and lower (/ =3, 4, 5) panel. Level of outermost grid is denoted in the 
upper left corner. The elapsed time tp and central density pc on the midplane are denoted 
in the top of each panel. The velocity scale in the unit of the sound speed is denoted in the 
bottom of each panel. 
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Fig. 2. — Evolution of average specific angular momenta derived in the region of r < 0.5 th 
for models with different finest grid levels Imax- 
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Fig. 3. — The average specific angular momenta jV against the distance from the origin with 
models having different sink radii fsink- The vertical dotted line represents the Hill radius. 
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3-75 



Fig. 4. — Structure around the Hill sphere for model M04 on the midplane (left) and in 
three dimension shown in bird's-eye view {right). The gas streamlines {red lines), density 
structure {color), and velocity vectors {arrows) are plotted in each panel. The dotted circle 
in the left panel represents the Hill radius. The size of the domain is shown in each panel. 
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Fig. 5. — Structure around the protoplanet for model M04 with different grid levels: / = 3 
(left), 5 {middle), and 7 (right). Each upper panel is the structure in three dimensions shown 
in bird's-eye view (right), and each lower panel is the structure on the cross section in the 
y = plane. The gas streamlines (red lines), density structure (color), and velocity vectors 
(arrows) are plotted in each panel. The size of the domain is shown in each panel. 
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Fig. 6. — Evolution of accumulated masses (panels a and 6), angular momenta (panels c 
and (i), and the average specific angular momenta (panels e and f) in the region of f < 0.1 
(panels a, c, and e) and f < 0.05 (panels 6, and f) against the elapsed time for models 
M005, MOl, M02, M04 and M06. 
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Fig. 7. — The density {false color, and contours), and velocity vectors (arrows) are plotted 
on the cross section in the z = plane in each upper panel for models M02 (left), M04 
[middle), and M06 (right). The Jacobi energy (false color, and contours) are plotted in each 
lower panel for the same models as in each upper panel. Three levels of grid (I =4, 5, and 
6) are superimposed in each panel. The dotted-circle represents the Hill radius. 




Fig. 8.— 



Same as Fig. [71 but in the cross section in the y 



= plane. 
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Fig. 9. — The distribution of accumulated mass (upper panel) and absolute value of the 
angular momentum (lower panel) against the distance from the protoplanet for models with 
MOOS, MOl, M02, M04, and M06. The thin solid lines represent the initial distribution of 
the mass (upper panel) and angular momentum (lower panel). The circles represent the Hill 
radii. 
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Fig. 10. — Average specific angular momentum jV against the distance from the protoplanet 
for models with MOl, M02, M04, M06, Ml, and M3. The circle and crosses represent the 
Hill and Jacobi radii, respectively. 
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Fig. 11. — Specific angular momenta in the region of r < rn (+), r < rj (o), r < 1/2 rn 
(A), and r < 1/10 th (□) against tlie protoplanet mass. Blue and red lines are the fitting 
formulae for M > IMj (blue) and M < IMj (red), respectively. 
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Fig. 12. — The ratio of the azimuthal to the Kepler velocity v^/vk {color and contours) 
around the protoplanet for models M02 (left), M04 (middle), and M06 (right). Three levels 
of grids are superimposed in each panel. The outer and inner white-dashed-circles represent 
the Hill radius f = fn (outer) and the half of the Hill radius f = l/2fH (inner), respectively. 
The black circle represents the contour of = Cs, inside which gas rotates with supersonic 
velocity > Cg. 



Fig. 13. — Density distribution {color and contours) on the cross section in the z = (upper 
panels), y = (middle panels), and a; = (lower panels) plane with different outermost grid 
levels / = 3 (left), 5 (right) and 7 (right) for model Ml. Three grid levels are superimposed 
in each panel. Right panels g, h, and i are eight times enlargement of left panels a, b, and c. 
The dotted circles represent the Hill radii f = fn. 



